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ABSTRACT 



Recently ISaito fc Sironol (|201lf) proposed that large ice aggregates which drift in- 
wards in protoplanetary disks break up during sublimation, ejecting embedded silicate 
particles. This would lead to a concentration of small solid particles close to the snow- 
line. In view of this model we carried out laboratory experiments where we observed 
freely levitating ice aggregates sublimating. We find that frequent break up is indeed 
very common. Scaled to a 10 cm aggregate about 2 x 10'' small silicate aggregates 
might result. This supports the idea that sublimation of drifting ice aggregates might 
locally increase the density of small dust (silicate) particles which might more easily 
be swept up by larger dust aggregates or trigger gravitational instability. Either way 
this might locally boost the formation of planetesimals at the snowline. 
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1 INTRODUCTION 

In protoplanetary disks small pa rticles collide, stick to gether 
and form larger aggreg ates (e.g.. lBlum fc Wurmllioog V This 
is usually regarded as a first phase of planet formation. 
One process with potentially wide impact on planet for- 
mation is radial transport of material. In a pressure gra- 
dient supported disk the disk rotates somewhat slower 
than Keplerian. As solids are not supported by this pres- 
sure gradient they feel a residual gravity and drift inwards 
(|Weidenschilling| (l977l V In typical disks meter-size bodies 
drift the fastest moving about 1 AU in 100 y. Usually, this 
fast radial drift is regarded as hurdle on the way to planets, 
as it might lead to the loss of solids to the star, eventu- 
ally. However, infall of mass from further out also provides 
means for a local concentration of solids further inward. Two 
mechanisms of t his kind, related to each other, were re cently 
proposed by Sait o fc Sironol (|201ll 'l and lSironol l|201lh . 

The basic idea proposed is that large aggregates drift 
inwards from the colder regions of the disk where the domi- 
nant mass fraction of solids is ice. They sublimate and during 
this phase fall apart in several pieces. Only the silicate cores 
embedded within the aggregates survive and due to the frag- 
mentation process a larger number of small particles or ag- 
gregates gather at or close to the snowline. As the fragments 
are smaller in size compared to their progenitor, the drift ve- 
locity is much smaller and the particles concentrate around 
the snowline. As consequence this might boost the forma- 
tion of planetesimals. The density might rise to a level that 
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gravit ational instability leads to the formation of planetes- 
imals (|Wakita fc Sekival 120081 '). Also in a turbulent disk an 
enhanced dust to gas ratio might be beneficial for planetesi- 
mal formation (|johansen et al.ll2009l ). Last not least growth 
efficiencies in the hit-and-stick scenario of planetesimal for- 
mation increase as well if larger bodies have a reservoir of 
small particles to collect (Blum and Wurm 2008). Whatever 
the details of further evolution, the breakup of large ice ag- 
gregates into smaller units while drifting inwards and the 
resulting concentration of small grains might be an impor- 
tant step in planetesimal and planet format ion close to the 
snowl ine. Therefore, the models proposed bv lSaito fc Sironol 
(|201ll ) and lSironol (|201ll ) offer interesting perspectives. The 
models have some distinctive features. 

The model bv ISironol (|201ll ) considers sublimation un- 
der saturated conditions including sintering. In a numeri- 
cal simulation the author shows how water is redistributed 
in an ice aggregate to the largest grain within this aggre- 
gate, eventually. During this process connections are broken 
and intermediately the aggregate is split into several smaller 
aggregates which eventually all sublimate completely. This 
model gets support from laboratory work where sublimat- 
ing / sintering snow layers under saturated condi tions are 
observed collapsing under gravity (|Flin et al.ll2003l ). 

The model bv lSaito fc Sironol l|201ll ) is similar but con- 
siders only sublimation of ice aggregates (no sintering) in 
a water vapor free environment. Here, the basic setup is 
that silicate cores are mantled by ice and these core-mantle 
grains form larger aggregates. Reaching the snowline the ice 
sublimates and all cores are assumed to be ejected as in- 
dividual grains. This model is not supported by numerical 
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calculations as the model bv lSironol (|201ll ). The fragmenta- 
tion of the ice aggregates itself is not discussed but a basic 
assumption. In view of this we carried out laboratory ex- 
periments for free particles suspen ded against grav i ty wh ich 
very closely resemble the model bv lSaito fc Sironol (|201ll ') to 
study how many fragments would result from sublimation. 

The experiments are based o n a levitation technique 
which was only recently developed. iKeUing fc WurmI (|2009l ) 
showed that porous dust particles levitate over a hot sur- 
face at low (mbar) pressure. This is due to thermal creep 
in a temperature gradient which transports gas downwards 
through the pores of the particles. This flow leads to a slight 
pressure increase below the particle which is sufficient to 
compensate its gravity. The same effect is used here but 
at lower temperatures. Details are given below. With this 
technique it is possible to observe free, levitating ice aggre- 
gates while they sublimate and break up which gives a direct 
measure of the number of times a given aggregate breaks up 
before it is completely sublimated. The aggregates in our 
experiments move is a quasi 2d way where gravity and com- 
pression due to gravity are only minor factors. 
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2 LABORATORY EXPERIMENTS 

We levitate ice aggregates at low pressure (mbar) which 
allows a sublimation of freely moving particles where the 
breakups of particle parts get directly observable and the 
number of breakups per time can be counted. The levitation 
of the particles is based on a n overpressure between aggre- 
gates and underlying surface l|Kelling fc Wurmll2009l '). This 
overpressure is generated by thermal creep through the pores 
of the particle. The thermal creep is induced by a temper- 
ature gradient and this temperature gradient is externally 
generated by cooling the top of the experiment chamber with 
liquid nitrogen. The pressure increas e by thermal cr eep is 
also known as Knudsen Compressor (|Knudsenlll909l ). The 
compression depends on the ambient pressure and allows a 
particle lift in the mbar pressure range. A particle will lev- 
itate at a height where the gas flow to the sides decreases 
the pressure difference to exactly the level that the weight 
of the particle is balanced, (see Fig. [T]) 

The actual setup consists of a Peltier element (5cm x 
5cm), on which an aluminum plate is placed. The plate is 
slightly concave to confine the aggregates to a limited area. 
Typical measured temperatures on the aluminum plate sur- 
face are 260 K. On top of the setup a reservoir of liquid 
nitrogen at 77 K is placed. A sketch of the setup can be 
seen in Fig. [2] 

At the low pressure, the ice sublimates at a convinient 
rate (within a few minutes). As seen in Fig.[T](b) the levitat- 
ing aggregates are typically fractions of a mm in thickness 
which is small compared to the total area of ice coverage of 
about 10 mm. This reduces the problem effectively to a 2D 
problem and breakups can well be observed from a top view 
of the ice particles (Fig. |3l Fig[4]). To do so, the nitrogen 
reservoir has a central clearance (not indicated in Fig. [2| . 

In the experiments gravity is well balanced by the pres- 
sure difference. Only small forces remain which can lead 
to an interaction between different particles, i.e. if contacts 
break aggregate parts are usually drifting apart. Due to the 
flat and porous structure vertical breakups with a succeeding 



Figure 1. a) Levitation induced by thermal creep (Knudsen com- 
pressor). The air flows through the pores from the cold to the 
warm side building up a pressure (Kelling and Wurm 2009). b) 
Side view of a levitated ice aggregate on a Peltier element and its 
reflection. 
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Figure 2. Sketch of the setup used for ice aggregate levitation. 
Experiments are carried out at a pressure of few mbar. 



compaction (aggregate parts fall downwards) are expected 
to be rare though we cannot rule out such events. We only 
consider observable breakups here, where one part of an ag- 
gregate clearly separates from another part. In each exper- 
iment a movie is taken from the sublimating ice particles. 
On the dark background a separation can easily be seen by 
eye while scanning frame by frame through the movie af- 
terwards. The number of breakups and the times of each 
breakup are the results of a visual analysis of the movie 
without any further image analysis. 

We also count breakups where two aggregate parts sep- 
arate but are still connected at another location, i.e. bend 
and open at the breakup point. However, such events only 
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Figure 4. Example of an observed aggregate breakup. In this 
sequence it can be observed in (a) how the lower part of the left 
aggregate is still connected to the upper part. In (b) the lower 
part has displaced slightly to the right but still remains conected. 
In (c) the aggregate is clearly broken and the lower part separates 
from the upper one. Finally (d) the two parts will continue their 
paths separately. 

account for about 10% of all cases. This is in agreement with 
the formation of the aggregates by hit-and-stick growth dur- 
ing preparation, where every individual ice grain only has 
two neighbors. 

Particle Preparation: Ice particles were prepared under 
atmospheric conditions (1 bar, laboratory humidity). Indi- 
vidual ice grains condensate in the vicinity of the pipe's sur- 
face that connects the liquid nitrogen tank and the setup's 
liquid nitrogen reservoir. These ice grains collide, stick to- 




Figure 5. Example from the preparation process. Atmospheric 
water vapor condenses into ice grains (small particles) around a 
cold liquid nitrogen pipe (top black area). They collide and stick 
to aggregates which stick to the pipe. 

gether and form aggregates, which can be observed by mi- 
croscopy. Fig. [5] shows a snapshot of the preparation phase. 
We did not specify the morphology or structure of the indi- 
vidual grains in more detail. The aggregates are taken from 
the pipe's frozen surface and put into the setup using a small 
spoon cooled by liquid nitrogen. Also, the transfer of parti- 
cles into the experiment chamber might induce modifications 
e.g. by direct condensation of water vapor but essentially the 
ice particles observed in the experiments are grown aggre- 
gates. 

As aggregates the particles are an appropriate model 
for large particles in protoplanetary disks which are als o 
assumed to be formed by aggregation (|Blum fc Wurml2008l ). 

An example of the evaporation rate is seen in Fig. [S] 
which shows the evolution of the evaporating ice during an 
individual experiment as fraction of the image covered by 
ice. The sublimation can well be described as linear with 
time but with a pronounced change in slope at about 25 s. 
Within the scope of this paper we consider either one or two 
linear approximations for all sublimation curves. We then 
get constant sublimation rates 

dA 

a = — — = const ( 1 ) 

where a is —1.25 ± 0.01 mm^/s initially (< 25 s in Fig. (Sjl 
and —0.41 ± 0.01 mm^/s for the later times in the given 
example. This is a difference of about a factor 3. 

For the breakup rate (number of observed breakups per 
second) we observe similar splits in two parts as for the 
evaportation rate which are correlated in time. This is shown 
in Fig. [3 

To quantify the total number of breakups for a given 
size we measure the breakup rate and the total time it takes 
the aggregate to sublimate from the initial size. This time 
is determined by the sublimation rate a and the inital par- 
ticle size. As initial particle size we considered the two di- 
mensional radius of gyration of the particle system at the 
beginning of an experiment defined as 

with Ti being the coordinates of the i-th pixel covered by ice 
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Figure 6. Evaporation for a selected experiment. Shown is the Figure 9. The breakup rate for ice particles with different ini- 

fraction of the image covered by ice over time. A change in slope tial size. Upper and lower values are correlated to the difference 

is visible at 25 s. iri sublimation rate. Each data point corresponds to a different 

experimental run. Solid lines mark average rates. 
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Figure 7. Number of breakups over time. Also here a linear de- 
pendence can be observed which is split in time in two parts. 
The split is correlated to the split in evaporation rate of the same 
experiment (see Fig.O 

and Vs being the coordinates of the center of mass and Np 
is the total number of pixels. 

Fig. [8] shows the ice surface coverage A over size r. The 
data follow a power law with power d = 2.0 ± 0.2. This is 
consistent with a porous but otherwise uniform coverage of 
ice in 2 dimensions independent of size or d = 2. Taking 
d — 2 the ice covered area can be approximated hy A = br'^ 
with h = 2.20 ± 0.26. 
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Figure 8. Area covered by ice over the size of an aggregate at 
the beginning of an experiment. Each data point represents an 
individual experiment. 



In the cases where we have to consider two distinct 
breakup rates we relate the second breakup rate to the size 
at the time the breakup rate changes (Fig. [7]). If a change 
in sublimation rate is not obvious (mostly for the smaller 
aggregates) only one average breakup rate and one size are 
considered. 

Fig. m shows the breakup rate dependence on size. Up- 
per and lower values are correlated to the difference in sub- 
limation rate. The average breakup rate differs by about a 
factor of 3 which is the same factor as for the correlated 
difference in evaporation rate. 

Fig. ID suggests that the breakup rate is not dependent 
on the aggregate size. We would like to caution though that 
the breakup constancy is not a universal feature in detail. 
As described above, aggregates form during the preparation 
of the ice particles. Individual grains have a size distribu- 
tion between several up to a few tens of microns, though 
the latter might already be small aggregates (Fig. [5]). This 
size distribution directly translates to a distribution in con- 
tact areas or necks between neighboring grains. We cannot 
spatially resolve these contact areas. For aggregates with 
a distribution of neck sizes the breakups proceed continu- 
ously in time with some statistical variations. Besides this 
the neck size distribution determines the breakup rate in 
detail. It is the inability to quantify this distribution that 
the breakup rate of all expenments is approximated by one 
average value. We note though that average breakup rates 
in individual experiments with similar initial size might dif- 
fer by a factor of 3 (see fig. (Qj! More complex breakup rate 
dependencies on size or even for a given size would certainly 
be possible within the uncertainties of the data. However, 
these simplifications allow us to specify a total number of 
breakups for a given size and scale this to much larger ag- 
gregates in protoplanetary disks. 



3 DISCUSSION 

The experiments give direct evidence that sublimating ice 
aggregates frequently break up as has been proposed by 
Saito and Sirono (2011) and - including sintering - by Sirono 
(2011). This will generate a number of smaller aggregates 
which eventually reduce to aggregates of embedded silicate 
(dust) particles. This is an interesting novel idea on the road 
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to explain the formation of planetesimals. Our laboratory 
experiments support this showing that breakup is indeed a 
viable mechanism to increase the density of smaller parti- 
cles in the vicinity of the snowline. We consider the ice in 
our experiments as suitable analog to cluster-cluster aggre- 
gates in protoplanetary disks. It is likely that some initial 
difference in ice structure and morphology is responsible for 
the two different sublimation rates observed. For each of 
the two existing sublimation branches we assign an aver- 
age breakup rate, n. The total number of breakups, N, is 
therefore NT, where T is the total sublimation time. The 
given average surface sublimation rate (on each of the two 
branches) a = dA/dt gives a linear total time T with aggre- 
gate surface, T = A/ \a\. \n terms of size this translates to 
&T — br^ / \a\. The total number of breakups then is 



As n scales with a factor 3 in the experiments correlated 
to the same scaling factor for a the actual sublimation rate 
is not important and we combine all empirical factors in 
a single constant parameter x — 2 mm~^va eq. O It had 
to be expected that the total number of breakups does not 
depend on time or any time dependent breakup rates or 
total sublimation rates. So the result that the total number 
of breakups only depends on the size is consistent with these 
expectations. 

The laboratory situation is effectively a 2d scenario 
where the 2d projection covered by ice is proportional to 
the square of the size. Ice aggregates in 3d have the same 
projected area if they have a fractal dimension of 2. This 
corresponds to aggregates which grow via cluster-cluster ag- 
gregation (CCAs) where aggregates of the same size hit 
and rigidly stick to each other at the first contact point. 
This growth process initially occurs in protoplanetary disks 
(Blum and Wurm 2008). For fractal dimension 2 the pro- 
jected surface essentially also scales with dimension 2. That 
means that the size of protoplanetary ice aggregates corre- 
lates well with the size of aggregates in our experiments. 
Therefore our results can directly be applied to large aggre- 
gates in protoplanetary disks and eq. [3] gives the number of 
breakups for a CCA cluster of given 3d size r. 

Assuming the largest ice aggregate in the form of a CCA 
would be about 10 cm in a protoplanetary disk then the 
aggregate would fragment into ~ 2 x 10'' pieces, a 1 cm 
CCA would fragment into ~ 200 pieces, releasing a large 
number of small dust aggregates, eventually. More compact 
ice aggregates will certainly fragment as well but the number 
of total breakups might vary. 

The laboratory experiments are closely resembling drift- 
ing particles in protoplanetary disks. Certainly, in detail the 
total immber of breakups might differ from the laboratory 
value. Nevertheless, eq. |3]might be a suitable approximation 
from the experimental side. Clearly not all dust particles em- 
bedded in an ice aggregate will leave as single particles as 
the shrinking ice aggregate only forms a limited number of 
smaller ice aggregates. We assume (not the result of this 
study) that every ice fragment with silicate cores will even- 
tually be reduced to an aggregate of its silicate cores with 
the same morphology, i.e. a CCA ice aggregate results in a 
CCA silicate aggregate. 



The prin ci ple idea behind the proposals by 
ISaito fc Sironol \2Qll\ ) and ISironol l|201lh is that the 
released silicate particles drift slower. Actually, for CCAs 
the overall aggregate size does not change the gas-grain- 
coupling time (Blum and Wurm 2008). As the gas-grain 
coupling time determines the drift velocity such CCAs 
would always drift with the same speed. This is due to 
the fact that the gas-grain coupling time is proportional 
to M/A, where M is the aggregate mass and A is the 
aggregate cross section or surface. For CCAs consisting 
of the same individual grains M/A is constant. However, 
the individual grain size of the cores is - per definition - 
smaller than the size of the individual ice grains. Also the 
density of silicates is about twice or three times as high as 
the density of water ice. This decreases the M/A ratio by a 
certain factor. Therefore, these final dust aggregates couple 
better to the gas of the disk and their radial drift velocity 
is indeed reduced. 

For compact aggregates (not subject to this study) 
breakup might also occur. If final silicate aggregates are com- 
pact the mass scales with M ~ r'^ and M/A depends linear 
on the overall size of the aggregates. Any smaller fragment 
then drifts slower than the original large ice aggregate. 

Therefore, while the details depend on the initial ice ag- 
gregate morphology and the details of the break up including 
silicate cores the break up of ice aggregates due to sublima- 
tion at the snowline will occur and it almost inevitably leads 
to a density enhancement of silicates. 



4 CONCLUSION 

Saito and Sirono (2011) and Sirono (2011) proposed that 
ice aggregates drifting inwards and sublimating at the snow 
line leave behind a large number of smaller silicate or dust 
grains with much smaller radial drift velocity. This leads to 
a dust density enhancement close to the snowline and subse- 
quently an increased likelihood for planetesimal formation. 
Planetesimal formation might be due to more frequent stick- 
ing collisions and growth or due to gravitational instabilities. 
The proposal by Saito and Sirono (2011) was based on the 
mere assumption that ice aggregates break up upon subli- 
mation and release silicate particles. We verified these ideas 
of ice aggregate breakup upon sublimation in laboratory ex- 
periments. Indeed we find that ice aggregates break up in a 
large number of smaller aggregates. Scaled to protoplanetary 
disks a 10cm CCA aggregate would fragment into ~ 2 x 10^ 
pieces. We only observed the fragmentation of pure ice ag- 
gregates. We currently assume that all silicate particles that 
would be embedded in one ice fragment would end up as 
one aggregated particle after total sublimation, i.e. a 10cm 
CCA ice aggregate would result in ~ 2 x 10* silicate aggre- 
gates. If the presence of silicate cores would influence the 
sublimation process or what the morphology of the silicate 
aggregates looks like is subject to future studies. However, 
the observed breakup into smaller fragments supports the 
idea that sublimation of ice aggregates close to the snow- 
line might promote the local formation of planetesimals by 
providing a reservoir of silicate particles that drift slower. 
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